Spin state transition in LaCoOa by variational cluster approximation 
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The variational cluster approximation is applied to the calculation of thermodynamical quantities 
and single-particle spectra of LaCo03 . Trial self-energies and the numerical value of the Luttinger- 
Ward functional are obtained by exact diagonalization of a C0O6 cluster. The VCA correctly 
predicts LaCoOs as a paramagnetic insulator and a gradual and relatively smooth increase of the 
occupation of high-spin Co 3+ ions causes the temperature dependence of entropy and magnetic 
susceptibility. The single particle spectral function agrees well with experiment, the experimentally 
observed temperature dependence of photoelectron spectra is reproduced satisfactorily. Remaining 
discrepancies with experiment highlight the importance of spin orbit coupling and local lattice 
relaxation. 

PACS numbers: 72.80.Ga,71.27.+a,79.60.-i,74.25.Ha 



I. INTRODUCTION 

LaCoOs has received considerable attention over the 
years because it seems to undergo two electronic transi- 
tions or crossovers in the temperature range between 50 
and 600 Kelvin 0,0. The first crossover, usually referred 
to as the spin state transition, can be seen most clearly 
in the magnetic susceptibility 0, IE]- Below 50 Kelvin 
LaCoOa is nonmagnetic, x ~ 0, indicating that all Co 3+ 
ions are in the low spin (LS) 1 A\ g or t\ g state realized for 
d 6 in cubic symmetry with sufficiently large crystalline 
electric field (CEF). Then \ rises sharply which indicates 
the thermal excitation of states with nonzero spin and af- 
ter a maximum around 100 Kelvin decreases again. In in- 
elastic magnetic neutron scattering Q the low frequency 
magnetic scattering intensity near T shows a very similar 
temperature dependence as x- A pronounced anomaly is 
also observed in the coefficient of thermal expansion [1,0, 
the heat capacity shows only a weak anomaly at the spin 
state transition^. 

Abbate et ai.[{J found that the valence band photoemis- 
sion spectrum (PES) and Ols X-ray absorbtion spec- 
trum (XAS) show little or no change across the spin 
state transition. On the other hand, Haverkort at al. [10t| 
found a significant temperature dependence of the Co- 
1^2,3 XAS and the X-ray magnetic circular dichroism 
spectrum (XMCD) below 500 Kelvin. Thornton et al.\l)\ 
observed a temperature dependence of the Co K-edge 
prepeak between 140 Kelvin and 800 Kelvin and Medarde 
at aL[T2] found that this temperature dependence sets in 
at 50 Kelvin i.e. the onset of the spin state transition. 
The nature of the spinful excited state which is respon- 
sible for the spin state transition has been under debate 
for some time. While it was proposed originally that 
this is the high spin (HS) 5 T2 g (or t\ g e 2 g ) excited state of 
the Co 3+ ion0, Korotin et aL[l3[ concluded from their 
LDA+U calculation that this state rather is the interme- 
diate spin (IS) 3 Ti g (or iig e g) state. Recently, however, 
experimental evidence has accumulated [1 fj 1 14. 1 1 5l 1 1 q. 1 1 7] 
that it is really the HS state which is populated. This 
leads to a certain puzzle in that a model calculation with 



an A% g ground state and a 5 T2 g excited state with fixed 
activation energy A = E( 5 T2 g ) — EQ-A\ g ) cannot repro- 
duce the experimental curve. If A is adjusted so as 
to reproduce the temperature where \ starts to deviate 
from zero - « 50 Kelvin - the resulting maximum value 
of x(T) near 100 Kelvin exceeds the experimental value 
by a factor of « 10. The fit is much better asuming an 
IS excited state, which has led some authors[j|0 to con- 
clude that an IS state is responsible for the transition. 
On the other hand, Haverkort et al. pointed out [10] that 
their XAS and XMCD spectra can only be explained by 
admixture of a 5 T2 g excited state and concluded that the 
activation energy E( 5 T2 g ) — E(}Ai g ) is temperature de- 
pendent, rising from 20 meV at 50 K to 80 meV at 700 
K[10( and leading to a much slower increase of the pop- 
ulation of HS ions with temperature. A somewhat puz- 
zling feature of this scenario is that in a situation where 
the ground state is LS an increase of the activation en- 
ergy implies an increase of the CEF splitting lODq with 
temperature. The increase [l8[ of the Co-0 bond length - 
the most important parameter determining lODq - with 
temperature, however, would result in exactly the oppo- 
site behaviour. The trend thus cannot be explained in a 
single- ion-pict ur e but is a kind of 'band effect'. Kyomen 
et a/.[H,[20|, who deduced a very similar temperature de- 
pendence of the activation energy to reconcile magnetic 
susceptibility and specific heat data, invoked a negative 
energy of mixing between low spin and high spin ions. 
Another reasonpLOj for the smaller-than-expected value 
of x observed in experiment is spin-orbit coupling which 
splits the 5 T2 g multiplet into three levels spanning an 
energy of w 75 mel/|16||. which is large compared to 
the temperature where the spin state transition occurs. 
The lowest of these spin orbit-split states is threefold de- 
generate. This low energy triplet - itself slightly split 
by the trigonal distortion due to the orthorhombic crys- 
tal structure of LaCoOs - can be identified in elec- 
tron spin resonance (ESR) [13, EH and inelastic neutron 
scattering [H, [2l| experiments, whereby its g-factor of 
~ 3 — 3.5 is clear proof that it originates from a spin- 
orbit-split HS state rather than from an IS state. 
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Haverkort et a/[lj| also found that to fit their XAS and 
XMCD spectra with cluster calculations they had to use 
a larger lODq for the A\ g state than for the 5 T2 9 state. 
This hints at a participation of lattice degrees of free- 
dom in that oxygen octahedra around HS Co 3+ expand 
slightly so as to accomodate the somewhat larger radius 
of the HS ion. The emerging picture thus is a disordered 
mixture of LS and HS states, whereby the lattice partici- 
patesby an expansion of the CoOg octahedra around HS 
sites [lOj, which would immediately explain the anomaly 
of the coefficient of thermal expansion 0, 0]- As pointed 
out by Berggold et aL[23| this idea also nicely explains an 
anomaly in the thermal conductivity n of LaCo03. At 
low temperatures the dominant contribution to K comes 
from phonons and the expanded 06 octahedra around 
HS Co constitute randomly distributed lattice imperfec- 
tions which reduce the mean free path of the phonons. 
This leads to a decrease of k at the onset of the spin state 
transition and a minimum slightly below 200 Kelvin. 
The idea of expanded 06 octahedra around HS Co ions 
may resolve yet another puzzle, namely the result from 
inelastic magnetic neutron scattering^, [2l| that low 
energy spin correlations in LaCo03 are ferromagnetic 
rather than antiferromagnetic. This is surprising in that 
the HS state has the configuration t\ g e 2 g so that the 
Goodenough-Kanamori rules would predict strong anti- 
ferromagnetic exchange interaction between two HS Co 
ions on nearest neighbors. The expansion of the oxy- 
gen octahedra around HS ions, however, would make the 
formation of HS states on nearest neighbor Co ions en- 
ergetically unfavourable, so that this antiferromagnetic 
nearest-neighbor exchange may never have the chance to 
act. The ferromagnetic spin correlations could then be 
due to 'semiconductor version' of the double exchange 
mechanism 22] . 

Strong experimental evidence against any appreciable oc- 
cupation of IS states is also provided by the EXAFS re- 
sults of Sundaram et al.\ll\. These authors ruled out 
the existence of inequivalent Co-0 bonds which would 
be practically inevitable in the presence of IS states be- 
cause the single electron in the two e g orbitals would 
make these strongly Jahn- Teller active. 
The second crossover in LaCo03 is frequently referred 
to as a metal-insulator-transition. It can be seen most 
clearlyin the specific heat where the raw data of St0len 
et aZ.jg] show a sharp 'spike' at 530 Kelvin even be- 
fore subtraction of the phonon background. Surprisingly 
for a metal-insulator-transition the electrical conductiv- 
ity a does not seem to show any noticeable anomaly 
at this temperature. Thornton et aL[2J] found that 
at low temperatures the electrical conductivity a shows 
a semiconductor-like increase with temperature which 
can be fitted well by assuming an activation energy of 
A = 0.53eU between 380 K and 520 K. There is a broad 
plateau between 600 K and 800 K and only above 800 
Kelvin a decreases with temperature as in a metal. Bhide 
et al. [H fitted the temperature dependence of a with an 
activation energy between 0.1 eV - 0.2eV for tempera- 



tures below 400 K. They found a plateau between 650 
K and 1000 K and a decrease with temperature only 
above 1200 K. Thornton et al^k, inferred a 'high-order 
semiconductor-to-metal transition' between 385 Kelvin 
and 570 Kelvin from a study of inflexion points in the a 
versus T plot. 

The magnetic susceptibility x nas a shallow maxi- 
mum near 600 Kelvin[l H" whereas magnetic neutron 
scattering @ does not show a pronounced signature of the 
transition. Abbate et aZ.Q found a significant change 
in the Ols XAS spectra between 100 Kelvin and 570 
Kelvin but the data of Thornton et al. (XAS at the 
Co K-cdgc) show a similar change as the Ols XAS as 
the temperature changes from 140 to 300 Kelvin so this 
change is not necessarily related to the metal-insulator- 
transition. Tokura et al. observed the filling of a gap- 
like structure in the optical conductivity <j(ui), Richter 
et al. did not observe a Fermi edge at temperatures 
above the crossover[27] in their photoemission spectra. 
The evidence for a true metal-insulator-tansition thus is 
not really compelling and in fact St0len et aL§| consid- 
ered an entirely different scenario, where the splitting of 
the 5 T2 g state by spin orbit coupling plays a central rule. 
In this scenario, the low temperature crossover is due to 
the thermal excitation of the low-energy triplet, whereas 
the 'metal-insulator-transition' corresponds the popula- 
tion of the remaining components. 

It has been argued that the crystal structure of LaCo03 
may play a role as well. In the LDA+U calculations 
of Korotin et al. [131 ] the structural change with increas- 
ing temperature is sufficient to induce a phase transi- 
tion between magnetic and nonmagnetic ground states. 
Quite generally density functional calculations show a 
strong sensitivity of ground state properties to structural 
parameters pjol l3l| . An experimental result result which 
directly shows the importance of lattice degrees of the 
lattice structure for the magnetic state of the Co 3+ ion 
is the ferromagnetism observed recently in LaCoOs thin 
films under tensile strain[28| Unlike bulk LaCoOs these 
thin films have temperature independent photoelectron 
spectra[29|. 

LaCoOs clearly is a difficult problem for any kind of elec- 
tronic structure calculation and has been studied by vari- 
ous methods during the last years: standard density func- 
tional thoryflcl, LDA+U or GGA+U[H EH, M, M, El 
and dynamical mean- field theory [35]. As already men- 
tioned LDA band structure calculations incorrectly pre- 
dict the material to be a metal in the paramagnetic state 
for both, the ideal perovskite structure and the true or- 
thorombic structure [30]. Combined photoemission and 
bremsstrahlung isochromat spectroscopy (BIS) data[36| 
indicate a gap in the electronic structure although its 
precise magnitude is difficult to pin down because the 
BIS spectrum shows a slow and almost linear increase of 
intensity with increasing energy. Together with the satel- 
lite structures observed in valence band photoemission^ 
this indicates the importance of electronic correlations 
and suggests that at low temperature the material is 
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actually a correlated insulator. From the above discus- 
sion it is moreover clear that a realistic description of 
the temperature dependence of the photoelectron spec- 
tra and magnetic susceptibility requires a correct descrip- 
tion of the multiplet structure of the Co 3+ ion, and its 
interplay with the crystalline electric field. On the other 
hand the relatively small gap indicates that covalency is 
strong so that band effects obviously are important as 
well. LaCo03 therefore appears as an interesting test 
case for the variational cluster approximation proposed 
by Potthoff [37]. This method generates trials self ener- 
gies in a finite cluster - an octahedral CoC>6 cluster in the 
present implementation - so that the interplay between 
multiplet structure and crystal field splitting can be eas- 
ily included. Being based on exact diagonalization rather 
than Quantum Monte Carlo the VCA can access low tem- 
peratures as necessary for the case of LaCo03. On the 
other hand, the present implementation is based on an 
LCAO-fit to the band structure whose necessarily limited 
accuracy makes it hard to quantitatively include the ef- 
fects of changes of the lattice. Therefore all calculations 
were carried out for a rigid lattice, which for simplic- 
ity was chosen to be the s.c. ideal Perovskite structure. 
Bearing in mind the scenario inferred by Haverkort et al. 
- an inhomogeneous lattice distortion with C0O6 octa- 
hedra expanding or contracting locally in response to the 
spin state of the Co-ion - it is quite obvious that a quanti- 
tative agreement with experiment cannot be expected for 
any calculation for a rigid lattice. A quantitative discus- 
sion of the temperature dependence of x moreover would 
require to include spin-orbit coupling which was omitted 
in the present study to simplify the calculations. Bearing 
this in mind we may expect the present calculation, with 
a rigid lattice and no spin orbit coupling will reach at 
best qualitative agreement with experiment. As will be 
shown below, however, this goal is indeed achieved. 



II. VARIATIONAL CLUSTER 
APPROXIMATION 

The quantity which is subject to variation in the varia- 
tional cluster approximation (VCA) is the electronic self- 
energy More precisely the VCA seeks for the best 
approximation to the self-energy S(w) of a lattice system 
amongst the subset of self-energies which can be repre- 
sented as exact self-energies of a given finite cluster. The 
VCA is based on an expression for the grand potential 
of an interacting many-Fermion system derived by Lut- 
tinger and Ward(38j. In a multi-band system where the 
Green's function G(k, uj), the noninteracting kinetic en- 
ergy t(k) and the self-energy E(k, to) for given energy 
lj and momentum k are matrices of dimension 2n x 2n, 
with n the number of orbitals in the unit cell, it reads [3^| 

= ~ E e ""° +ln det (-G-^k,^) ) + F[E](l) 



where u v = (2v + l)ir//3 with j3 the inverse temperature 
are the Fermionic Matsubara frequencies, 

G-^k,^) = iu + h -t(k) - S(k,w). (2) 

with [i the chemical potential and F[E] is the Legendre 
transform of the Luttinger-Ward functional 3>[G]. The 
latter is defined [38] as the sum of all closed linked skele- 
ton diagrams with the non-interacting Green's functions 
replaced by the full Green's functions. A nonperturbative 
derivation of a functional with the same properties as $ 
has been given by Potthoff [4(| . Luttinger and Ward have 
shown that Q is stationary with with respect to variations 
of E: 



<9Ey(k,w„) 

but the crucial obstacle in exploiting this stationarity 
property in a variational scheme for E is the evaluation 
of the functional F[H] for a given 'trial E'. Potthoff's 
solution[37} makes use of the fact that -F[E] has no ex- 
plicit dependence on the single-particle terms of H and 
therefore is the same functional of E for any two sys- 
tems with the same interaction part of the Hamiltonian. 
In the following only the Coulomb interaction within the 
Co3c?-shell is taken into account - which is a reasonable 
approximation. Under this assumption F[S\ then is the 
same functional for the true perovskite lattice and for 
an array of identical but disconnected octahedral CoOg 
clusters. For given value of fj, and /3 one can therefore con- 
struct trial self-energies S(w) by exact diagonalization of 
a single Co06 cluster and at the same time obtain the 
exact numerical value of F[S\ by simply reverting the 
expression ([I]). Here the kinetic energy t of the Co06 
cluster has to be used in the Dyson equation ©. Next, 
the pair (E, F[E]) can be used in (JTJ) for the lattice sys- 
tem - which simply amounts to replacing t by the kinetic 
energy of the lattice, t(k), in ([2]) and performing the k- 
summation - to obtain an approximation for the grand 
potential of the lattice. The variation of E(w) then is 
performed by varying the single-electron parameters 
- such as hybridization integrals or site-energies - of the 
Co06 cluster. The condition © thus is replaced by the 
set of conditions 

I- 

Potthoff has introduced the name 'reference system' for 
the finite cluster used to construct trial self-energies and 
computing the Luttinger-Ward functional. In the present 
application - described in detail in Refs. [12] and [43[ - 
this is an octahedral CoOg cluster. Since it is known 
that exact diagonalization of clusters comprising a sin- 
gle transition metal ion and its nearest neighbor oxygens 
gives excellent results for the k-integrated photoelectron 
spectra of many transition metal oxides [3, |4f| [4(| one 
may expect that the use of such a cluster as the reference 
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system is a reasonable choice. However, different imple- 
mentations of the VCA have used quite different refer- 
ence systems. After being proposed by Potthoff[33] (an 
excellent review covering many technical details has been 
given by Senechal[41() the VCA has been applied succes- 
fully to one- and two-band Hubbard models [13, EH, to 
simplified models for to Fe pnictides [4^l a nd to transition 
metal oxides with orbital degeneracy [5fJ. 
To obtain the single-electron Hamiltonian t(k) an LDA 
band structure calculation for LaCo03 was performed 
using the Stuttgart LMTO-package. Thereby the ideal 
cubic perovskite structure with a Co-0 bond length 
of 1.91 A was assumed. The density of states is con- 
sistent with previous results and actually quite sim- 
ilar to that obtained for the correct rhombohedral 
structure [131 . Next, an LC AO-fit was performed to ob- 
tain a multi-orbital tight-binding parameterization of the 
single-electron Hamiltonian t(k). The LCAO basis com- 
prises 02s and 02p orbitals at —12.834 eV and eV, 
Co 4s and 3d orbitals at 19.436 eV and 1.731 eV and La 
5p and 5d orbitals at -9.264 eV and 10.436 eV. For the 
Co3d orbitals an additional lQDq = 0.848 eV is obtained 
from the fit. All orbitals except 02p and Co3d only help 
to 'polish' certain portions of the band structure, but to 
obtain a good fit they have to be included. The energies 
of these 'auxilliary orbitals' were not subject to the fit 
only the respective two center integrals. These are listed 
in Table |TJ In general they refer only to nearest neighbor 
bonds, but hybridization between second nearest neigh- 
bor oxygen has also been used. Figure [T] compares the 



(spa) 
(ppa) 
(ppn) 
(sda) 
ipda) 
(pdn) 



Co-O O-O O-O La-0 



-1.504 
0.000 
0.000 
-1.201 
1.776 
-0.975 



0.000 0.000 

0.930 0.132 

■0.112 0.000 

0.000 0.000 

0.000 0.000 

0.000 0.000 



0.000 
1.488 
-0.289 
0.000 
-0.879 
0.296 



TABLE I: Two center integrals (in eV) obtained by a LCAO 
fit to paramagnetic LDA band structure of LaCoOa 



actual LDA band structure and the LCAO-fit - while 
the fit is not really excellent the overall band structure 
is reproduced reasonably well. The bands at the top of 
the figure which are absent in the LCAO band structure 
originate from La5c? orbitals and would require to adjust 
parameters such as La (dda) - here they are put to zero 
because these bands are not really interesting. In the in- 
tervalls —10 eV — > —7 eV and —4 eV — ► eV there are 
bonding/antibonding bands of mixed Co3e?/O2p charac- 
ter, in the intermediate energy range there are essentially 
nonbonding 02p bands with very little Co3d admixture. 
The deviation from LDA seems large at the R point, 
but it should be noted that the band top at M is also 
somewhat lower in the LCAO-fit (the two band struc- 
tures were alligned at T). The total bandwith between 
the maximum at M and the minimum at R is 9.48 eV 




FIG. 1: (Color online) LDA band structure (left) and LCAO- 
fit (right) for LaCoOs in the ideal Perovskite structure. 



for the LDA and 9.83 eV for the LCAO-fit, i.e. the dif- 
ference is 4% which seems tolerable. 

The Coulomb interaction between Co 3d electrons is 
described by standard atomic multiplet theory [Fl|, 
More precisely the Coulomb interaction within the d-shell 
can be written as 



E 



(KlK 2 \g\K4K 3 ) 4i4 2 <i «3 < W 



(5) 



Here we have suppressed the site label i and k = (m, a) 
where to G { — 2,— 1,...2} denotes the z component of or- 
bital angular momentum. The Coulomb matrix elements 
(K1K2I3IK4K3) are obtained by a multipole expansion of 
the Coulomb interaction term l/|r — r| and involve Gaunt 
coefficients from the angular integrations and the three 
Slater-Condon parameters F°, F 2 and F 4 from the ra- 
dial integrations. The upper index thereby refers to the 
multipole order of the interaction and in a ci-shell is lim- 
ited to < 4 by the triangle condition. The somewhat 
lengthy complete expression for the matrix elements is 
given e.g. in equations (13-18)-(13-25) of the textbook 
by Slater [5l|. F 2 and F 4 which describe higher multipole 
interactions can be calculated from atomic Hartree-Fock 
wave functions but F° is reduced substantially from its 
atomic value by solid state screening and is treated as 
an adjustable parameter. In the present work the val- 
ues F° = 8.376 eV, F 2 = 10.64 eV" and F 4 = 6.804 eV 
or, alternatively, the 3 Racah-parameters, A — 7.62 eV, 
B = 0.14 eV, C = 0.54 eV were used. For the lowest 
mulpiplets the full theory as described by ([5]) can be re- 
duced to a parameterization in terms of a Hubbard U and 
Hund's rule </, which parameters then can be expressed in 
terms of the Slater-Condon parameters [Hj]. As discussed 
in Ref. [43J we also need to specify the 'bare' (i-level en- 
ergy e^. While the LCAO fit does give an energy td for 
the Co3d level, this contains a large contribution from the 
intra-d-shell Coulomb interaction. Since we want to de- 
scribe this Coulomb interaction by adding the Coulomb 
Hamiltonian to the LCAO-like single-particle Hamil- 
tonian we have to correct for this to avoid double count- 
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ing. For example, Kunes et aL[54( have estimated this 
double counting correction as e* d = — 9Un d where 
is the average electron number/d-orbital. With a U of 
order 10 eV this correction obviously is large. Since first- 
principle calculations of screened interaction parameters 
are a subtle issue, however, e* d was considered as an ad- 
justable parameter in the present work and set to be 
e* d = —46.4 eV. F° and e* d together essentially determine 
the magnitude of the insulating gap and the distance of 
the 'satellite' in the photoemission spectrum from the 
valence band top. With these values of the Racah pa- 
rameters and e* d and using Table III of Ref. [43[ we ob- 
tain the energy differences E(d n+1 ) + E(d n - 1 )-2E(d n )= 
A-8B = 6.5 eV E(d n+1 L) - E(d n )= 1.98 eV which are 
frequently referred to as the Hubbard U and the charge 
transfer energy A. Korotin et q/.[l3| obtained the value 
U = 7.5 eV by density functional calculations. 
Finally, we discuss the CEF splitting lODq. It is obvious, 
that the CEF is a crucial parameter for LaCo03 because 
it determines - amongst others - the relative energy of 
the A\g and 5 T2 g state of the Co(3+) ion. One can not 
expect that the LDA calculation and the LCAO-fit will 
produce a sufficiently accurate estimate so as to repro- 
duce energy scales of the order 100 Kelvin. Therefore 
lODq was also treated as an adjustable parameter and 
to get agreement with experiment the value 10Dq=0.72 
eV was chosen, which still is rather close to the value of 
0.848 eV obtained from the LCAO fit. 
Next we briefly comment on the technical problem of 
finding a stationary point of 11 in a multi-dimensional 
parameter space. As a first step, all but one parameter 
Ao are kept fixed and Ao is varied until a value where 
dfl/dXo = is found. This means we are now on a sur- 
face in parameter space - which we call the '(Ao)-surface' 
- where dfl/dX = 0. It is advantageous to always choose 
Ao to be the center of gravity of all orbital energies in the 
reference system because Aichhorn et al. have shown |48|. 
that optimization of this parameter leads to a thermo- 
dynamically consistent particle number. Then a sec- 
ond parameter Ai is chosen and varied. In each step 
Ao is recalculated to maintain dfl/dXo = i.e. we walk 
along the (Ao)-surface in parameter space while varying 
Ai. The recalculation of Ao can be done by means of 
the Newton-method, thereby using the solution for the 
preceding value of Ai as initial guess for the next one. 
Variation of Ai is continued until a value is found where 
dtt/dXi — 0. This means we have now found a point of 
the '(Ao, Ai)-surface ' in parameter space which is defined 
by dfl/dXo = and simultaneously dfl/dXi = 0. Next, 
we choose a third parameter, A2 and vary this again, 
walking along the (Ao, Ai)-surface - the recalculation of 
Ao and Ai in each step is again done by the Newton 
method - until we find a point where 8^1/8X2 = and so 
on. This method has the advantage that it is in principle 
guaranteed to find a stationary point. Moreover by do- 
ing a wider scan of one parameter one can find different 
branches of the A-surfaces which correspond to different 
stationary points. 



Finally we comment on the choice of the parameters 
to be optimized. Using the notation of Ref. [43| the 
4 parameters to, e\, £2 and V(e g ) were varied. The 
values for the remaining parameters, 63 = 1.4 eV and 
V(t2g) — 2(pd7r)[44| were kept fixed. It was checked that 
optimization of more than 4 parameters led to negligible 
change of fl and very small changes in the single particle 
spectral function. The reason for this 'saturation' of f2 
is the existence of 'nearly stationar y' l ines in parameter 
space as discussed in detail in Ref. [43j . 



III. RESULTS 

A search in parameter space for stationary points (SP) 
of Q revealed that for most temperatures there are actu- 
ally three different SP corresponding to a d-shell occupa- 
tion of « 6. At low temperature the first one corresponds 
to the reference system being in the pure Ai ff state, the 
second SP corresponds to the reference system being in 
a thermal mixture of an Ai g ground state and a 5 T2g 
state at slightly higher energy, whereas the third SP cor- 
responds to the reference system being in the pure 5 T2g 
state. For 'reasonable' parameters the third SP - corre- 
sponding to the HS state of the reference system - has 
an £1 that is substantially higher than that of the first 
two SP, whence this SP will never be realized. At low 
temperatures, on the other hand, the Ai g -like SP and 
the 'mixed' SP have very similar values of f2 and for suit- 
able choice of 10 Dq in the physical system a crossover 
can be seen. To illustrate this, Figure [2] shows O for the 
Ai g -like and mixed SP as a function of temperature for 
lODq = 0.72 eV (the value of lODq obtained from the 
fit to the LMTO band structure was 0.85 eV). Accord- 
ingly this value of 10-Dg was kept fixed for the rest of 
the calculation. Much unlike the cases of NiO, CoO and 
MnO[43j], the parameters of the stationary points have a 
strong temperature dependence in the case of LaCoC>3. 
Obviously this reflects the subtle change with tempera- 
ture of the electronic structure. 

One can recognize in Figure [21 that there is a crossing 
of the two f2(T) curves at w 50 Kelvin, and the finite 
difference in slope implies a first order phase transition. 
This way of describing the spin state transition is prob- 
ably an artefact of the mean-field-like description in the 
framework of the VCA. The latent heat for the transition 
would be TAS = 27.8 J/mol. It should also be noted, 
that the VL{T) curve of the Ai ff -like SP has an unphysi- 
cal upward curvature above « 50 Kelvin. The lower SP, 
however, does indeed have the correct downward curva- 
ture. 

More interesting is the entropy S(T) because this can be 
compared directly to experiment. Figure shows S(T) 
from the VCA and the experimental electronic entropy as 
extracted by St0len et al. [8fl from their specific heat mea- 
surements. At low temperatures the agreement is not so 
bad but it is immediately obvious that the crossover at 
530 Kelvin which is very pronounced in the entropy is 
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FIG. 2: (Color online) Left: Grand canonical potential Q for 
the two low energy stationary points as a function of temper- 
ature. The value of lODq = 0.72 eV in the lattice system. 
Right: Entropy for the 'mixed' SP as a function of tempera- 
ture compared to the experimental electronic entropy found 
by St0len et aZ.||. 
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FIG. 3: (Color online) Magnetic susceptibility of LaCoOa 
from Refs. 4, 5] and spin susceptibility obtained by the VCA. 
Note that the VCA result is multiplied by a factor of 4. The 
spin susceptibility in the 'Ai 3 -like SP, which is realized below 
50 Kelvin, is essentially zero. 



not reproduced at all by the present calculation. This 
will be discussed below. 

It should also be mentioned that the IS (or 3 Ti g ) state 
has negligible weight in the reference system even at 
the highest temperatures studied. The reason is sim- 
ply the fact - already noted by Haverkort et al. 10] - 
that the IS state never comes even close to the ground 
state of the octahedral CoC>6 cluster. In that sense, 
the VCA complies with 'LS-HS scenario' supported by 
experiment. 0, EI EE EE E3 

Next, Figure [3] compares the spin susceptibility from 
the VCA calculation to the experimental magnetic sus- 
ceptibility. While the overall behaviour is similar, the 
temperature where x has its maximum does not agree 
with experiment and, more importantly, the calculated 
values are a factor » 4 too small. It is likely that the 
reason is the rigid lattice used in the present calculation: 
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FIG. 4: (Color online) Occupation of the 5 T2 ff state in 
LaCoOa as inferred by Haverkort et aZ. [T3] from their XAS 
spectra and by Kyomen et al. [2(j from a fit to the suscepti- 
bility and specific heat compared to the 5 T2 9 occupation in 
the reference system for the VCA calculation. Also shown 
is the 5 T2 g occupation obtained with a fixed activation en- 
ergy of 250 Kelvin, which reproduces the onset the spin state 
transition^. 



in the actual material the expansion of the Oq octahedra 
around a Co-ion in the high-spin state probably prevents 
HS ions from occupying nearest neighbors, so that the 
strong antiferromagnetic superexchange cannot act. In 
the VCA calculation this effect is absent, whence the an- 
tiferromagnetic nearest-neighbor-exchange probably re- 
duces the ferromagnetic spin polarization induced by the 
magnetic field. 

One may ask for the fraction of Co-ions being in the HS 
state. It should be noted, that the VCA does not give 
that number for the physical system. The exact diago- 
nalization of the reference system does give the occupa- 
tion numbers of the different eigenstates of the reference 
system, but there is no justification for identifying these 
with the occupation numbers in the physical system. On 
the other hand, if the optimal self-energy for the lattice is 
realized in a cluster where the the HS state has a certain 
weight a it is reasonably plausible that the occupation of 
the HS state in the physical system will not be differ com- 
pletely from a. Thus, we may consider the occupation of 
the HS state in the reference system as a plausible esti- 
mate for the true HS occupation in the physical system. 
Figure U compares this number to experimental values. 
Most importantly the increase of the HS occupation with 
temperature is much weaker than for a system with fixed 
activation energy. The estimates of Haverkort et al and 
Kyomen et al. are reasonably close and also the VCA 
gives a roughly correct description although it obviously 
underestimates the HS population. This is another rea- 
son why the susceptibility computed by the VCA is too 
small. It is interesting, however, that the VCA gives the 
temperature dependence of the HS occupation at least 
qualitatively correct as it was carried out with a rigid 
lattice and therefore includes 'band effects' but no ef- 
fects of the local lattice relaxation. 
Next we consider the occupation numbers of the various 
orbitals. These can be obtained for the lattice system 
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FIG. 5: (Color online) Change of occupation numbers with 
temperature for the 'mixed' SP. The values at 50 kelvin have 
been subtracted off to make changes more clearly visible. Also 
shown is the change of occupation numbers plotted versus HS 
occupation in the reference system. 



in the standard way by integrating the spectral function 
of the lattice system up to the chemical potential. At 
50 Kelvin, the occupation numbers/spin direction and 
atom for the various orbitals are n(Co t2 g ) = 2.956, 
n(Co e g ) = 0.473, n(02p) = 2.823. It is immediately 
obvious from these numbers that there is considerable 
charge transfer from Oxygen to the e g orbitals of Cobalt. 
As HS Co is admixed with increasing temperature, the 
occupation numbers change, as can be seen in Figure [5] 
The figure also shows the change of occupation numbers 
plotted versus the HS occupation in the reference sys- 
tem, i.e. the quantity which is shown as 'VCA' in Figure 
2J The very accurate linear dependence is an indication 
that the HS occupation in the reference system is indeed 
a very good estimate for the HS occupation in the lat- 
tice system. The changes are as expected, with n(Co e g ) 
increasing at the expense of n(Co t2 g ) and a slight net 
charge transfer from Co to O. The nearly equal and op- 
posite change of n(Co t%g) and n(Co e g )is expected if 
HS ^2g e g are admixed to LS t\ g . Since the e g orbitals 
hybridize with O by the stronger a bonds and the ti g by 
the weaker tt bonds, admixture of HS states will decrease 
the degree of covalency hence the slight charge transfer 
back to Oxygen. More precisely, each additional HS Co 
transfers 0.073 electrons to the 02p bands. 
Next, we consider the single-particle Green's function cal- 
culated with the optimal self-energy. Figure [5] shows the 
k-integrated spectral function near the Fermi energy as 
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FIG. 6: (Color online) Combined PES and BIS spectra (from 
Ref. [36l]) compared to the k-integrated single particle spectral 
function obtained from the VCA at 100 Kelvin. 5-peaks are 
replaced by Lorentzians with a width of 0.02 eV. 



well as the combined PES and BIS spectra by Chainani 
at al. [36J] . Most importantly, the VCA correctly describes 
LaCoOa as a paramagnetic insulator - there is a clear gap 
of « 1 eV in the spectrum. The BIS spectrum does not 
really show a clear edge but a gradual increase so it is 
hard to deduce a unique experimental gap value. The gap 
in the VCA spectrum of « 1 eV is larger than the gap 
values deduced from the temperature dependence of the 
conductivity [3, l2~ij which range from 0.1 eV to 0.53 eV 
(depending on temperature) or from the optical conduc- 
tivity, 0.1 eV[4|. On the other hand it should be noted 
that the theoretical gap value does not have much real 
significance anyway - it is largely determined by the ad- 
justable parameters U and A. 

More specific is the overall shape of the photoemission 
spectrum. Figure [7] compares the k-integrated spectral 
function with the experimental photoemission spectrum 
over a wider energy range. More precisely, we consider 
the 'On-off-difference', that means the difference of va- 
lence band photoemission spectra taken with photon en- 
ergies on (63.5 eV) and slightly off (60.0 eV) the Co3p — > 
3d threshold, a procedure which is known Q to emphasize 
the Co3d-derived features. Also shown is a photoemis- 
sion spectrum taken with a photon energy of 21.2 eV 
where, due to the larger photoionization cross section of 
the 02p orbital at this energy, mostly the oxygen de- 
rived states are visible. The VCA reproduces the main 
features quite accurately: the high intensity d-like peak 
at —1 eV, the smaller d-like peak at —5 eV and the broad 
'satellite' around —12 eV. Also, the 3 02p-like peaks are 
reproduced. The fact that the peak at the top of the va- 
lence band has predominantly Co3d character also comes 
out correctly. By and large the VCA gives a reasonable 
description of the electronic structure of LaCoOa, at least 
on coarse energy scales. 

An interesting feature of LaCoOa is the temperature 
dependence of its photoelectron spectra and the VCA re- 
produces these at least qualitatively. Figure [5] shows the 
k-integrated spectral function at two different tempera- 
tures, 80 Kelvin and 570 Kelvin. Also shown in the inset 
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FIG. 7: (Color online) Angle integrated valence band photoe- 
mission spectrum obtained from the VCA for the 'mixed' SP 
at 50 Kelvin compared to the 'on-off ' spectrum and a photoe- 
mission spectrum taken at hv = 21.2 eV (Experimental data 
taken from Saitho et aZjjJ, Lorentzian broadening of the VCA 
specta: 0.1 eV). 
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FIG. 8: (Color online) Angle integrated valence band photoe- 
misssion spectrum of LaCo03 obtained by VCA at 80 Kelvin 
(bottom) and 570 Kelvin (top) (Lorentzian broadening 0.1 
eV). The horizontal line is a guide to the eye and corresponds 
to the same intensity in both spectra. The inset shows exper- 
imental spectra at different temperatures (Ref. 0). 



are experimental angle integrated photoemission spectra 
by Abbate et al.. As one can see in experiment the promi- 
nent peak at the top of the valence band looses weight 
with increasing temperature, which is a manifestation of 
the increasing number of Co-ions in the high spin state. 
The VCA reproduces this effect qualitatively, but higher 
temperatures are needed to obtain a similar degree of 
spectral weight loss. As discussed above this is simply 
due to the fact that the VCA underestimates the HS oc- 
cupation. 

Next, we consider the unoccupied part of the spectrum 
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FIG. 9: (Color online) k-integrated single particle spectral 
function from VCA at different temperatures (Lorentzian 
broadening 0.1 eV) compared to Co K-edge XAS spectra of 
LaCoOa (taken from Ref. [ll|]). The XAS spectra have been 
shifted downward by » 7 keV. 



and compare to the Co K-edge spectra by Thornton et 
a$n\ . Figure [9] shows the spectral function at several 
temperatures as well as the 'prepeak' of the Co K-edge 
spectra - shifted in energy so as to match the electron 
addition spectrum. The VCA spectra show a peak at 
rs 0.7 EeV above the Fermi energy which increases in 
intensity as the temperature increases. The experimen- 
tal spectrum shows a similar change, i.e. the growing in 
intensity of a low energy peak. Again, since the VCA un- 
derestimates the HS population the growth of this peak 
is probably underestimated. It should also be noted that 
the Ols XAS spectra of Abbate et al.^ show quite a 
similar behaviour as the Co K-edge spectra, namely the 
growth in spectral weight of a low energy peak with in- 
creasing temperature. An apparent difference between 
the results of Thornton et al. and Abbate et al. is 
that the Co K-edge spectra seem to show a rather grad- 
ual change of the spectra with temperature, whereas the 
Ols XAS spectra show little change at low temperature. 
The VCA gives a very continuous and gradual change as 
would be expected from the rather smooth increase of 
the HS population. 

Lastly, we proceed to a comparison with other calcu- 
lations on LaCoOa. Several authors have calculated the 
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FIG. 10: (Color online) k-integrated single particle spectrum 
obtained from the VCA for the 'Ai g -like' SP at 10 Kelvin, 
Lorentzian broadening 0.1 eV). 

density of states (DOS) for the low spin - or nonmagnetic 
- state using LDA+U or GGA+U 03 M, S3, H 1| and 
it may be interesting to compare the VCA to these calcu- 
lations. Figure [TO] shows the k-integrated single particle 
spectrum at the lowest temperature studied, 10 Kelvin. 
Whereas all spectra shown so far corresponded to the 
'mixed SP' in FigureQ] this spectrum is calculated for the 
'pure Ai 9 ' SP. Despite this, one can see that the spectrum 
is nearly indistinguishable from the other spectra, which 
shows that the phase transition from the pure Ai g SP 
to the mixed SP in Figure [T] has practically no influence 
on the spectrum. With the exception of the Co3d-likc 
'satelite' at « —12 eV the spectrum is quite consistent 
with the GGA+U calculation of Pandey et al. [32| . Espe- 
cially the respective oxygen or Co character of the three 
prominent peaks agrees reasonably well and these agree 
in turn with the photon-energy dependence of the PES 
spectraQ. The DOS obtained by Hsu et al.^3^ shows 
three prominent peaks as well, but the characters do not 
match: there, the topmost peak has predominant oxygen 
character, whereas the lowermost peak has predominant 
Co character - this does not agree with experiment. 

IV. DISCUSSION 

To summarize one may say that the VCA gives a rea- 
sonably accurate description of some experimental results 
for LaCoOs. The insulating nature of the material is 
described correctly and the photoelectron spectra agree 
with experiment in quite some detail. The temperature 
dependence of the magnetic susceptibility and the photo- 
electron spectra is reproduced at least qualitatively. For 
all temperatures studied only LS and HS states have ap- 
preciable weight in the density matrix of the reference 
system, which means that the VCA is consistent with 
the LS-HS scenario supported by experiment. Thereby 
the population of the HS states increases quite smoothly 
with an onset at 50 Kelvin which would be consistent 
with experiment as well. 

The main deficiencies are the failure to reproduce the 



crossover seen at 530 Kelvin in the entropy and suscepti- 
bility, the too small value of the magnetic susceptibility 
and the too slow increase of the HS population with tem- 
perature, which makes the temperature dependence of all 
photoelectron spectra weaker than observed. It should be 
noted that changing the values of A and/or U so as to 
obtain e.g. a smaller insulating gap does not change this. 
The above deficiencies are very probably not related to 
an inappropriate choice of parameters. 
The first reason for deviations is probably the neglect of 
spin-orbit-coupling in the d-shell. This leads to a split- 
ting of the 5 T2 g state into a three-fold, a five-fold and 
a seven- fold degenerate state [13, EH, EU which span an 
energy of 75 meV. If one were to assume that the acti- 
vation energy in the reference system as obtained by the 
VCA corresponds to the center of gravity of these split 
states, the triplet would be appreciably lower and give a 
higher susceptibility and HS occupation at low temper- 
atures. Moreover, St0len et aZ.[|[ have considered a sce- 
nario, where the low temperature crossover corresponds 
to the population of the low energy triplet and the high 
temperature transition to the population of the remain- 
ing two components. If that were indeed the case, a cal- 
culation without spin orbit coupling could never repro- 
duce the high temperature transition. Since spin-orbit 
coupling is a single-particle term it can be included into 
the VCA without any problem. On the other hand the z- 
component of the spin is no longer a good quantum num- 
ber if spin-orbit-coupling is introduced, which increases 
the size of matrices to be diagonalized or inverted and 
spin-orbit coupling was neglected in the present study. 
A second reason is the neglect of - or rather: the im- 
possibility to treat - the local lattice relaxations i.e. the 
expansion of Oq octahedra around HS Co ions. This im- 
plies than a HS ion 'feels' a different environment and 
that the actual activation energy has to be modified by 
an elastic contribution. All of these effects are missed in 
a calculation with a rigid lattice like the present one. The 
local expansion of the 06 octahedra may also have con- 
siderable impact on the magnetic susceptibility. Namely 
based on the Goodenough-Kanamori rules HS ions on 
nearest neighbors would be expected to show strong anti- 
ferromagnetic exchange via the two half-filled e g orbitals. 
On the other hand, for HS ions on nearest neighbors 
it is clear that the respective local lattice relaxations - 
expansion of the Og octahedra around HS ions - would 
interfere with each other, so that HS occupation of near- 
est neighbors may be energetically unfavourable and the 
antifcrromagnetic superexchange may simply not have 
the chance to act. This could explain the experimental 
result @, [2lJ that low energy spin correlations are ferro- 
magnetic rather than antiferromagnetic as well as the 
surprising fact that thin films of LaCoOs under tensile 
order ferromagneticallv [28j . In a calculation with a rigid 
lattice this effect would be missed, so that the antiferro- 
magnetic superexchange would reduce the spin suscepti- 
bility. 

This shows that important physical effects had to be 
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neglected in the present calculation and a quantitative 
agreement with experiment could not be expected. Still 
there is quite good qualitative agreement which demon- 
strates the usefulness of the VCA to study correlated 
insulators. 

I would like to thank K. P. Bohncn, D. Fuchs, M. 
Haverkort, M. Potthoff and S. Schuppler for instructive 
discussions. 



Taking B small but finite and demanding that J^- 
we obtain 



£ 



d 2 n 

dXidX. 



d 2 tt 



We now differentiate with respect to B and set B = to 
obtain 



V. APPENDIX: CALCULATION OF THE 
MAGNETIC SUSCEPTIBILITIY 

The magnetic susceptibility can be obtained from x = 
— §W- Thereby the magnetic field B is an additional 
single-particle-like parameter in the physical system. The 
introduction of this parameter will change the stationary 
point, that means the parameters of the reference system 
become depedent on B. For the calculation of the deriva- 
tive, however, we do not need to solve the optimization 
with applied B-field. 

We denote by Xi the parameters of the reference system 
and by Xi the values at the stationary point for B = 0. 
Then we can write down the following expansion of Q for 
small B: 



n 



d 2 n 

dXidX 



d 2 n 

dXidB 



(A i -A i )(A J -A J ) 

B 2 d 2 Vt 



B(Xi - Xi) + 



2 dB 2 



(6) 



All derivatives in this equation can in principle be ob- 
tained numerically. Also it has been used that in a para- 
magnetic state 



on 

~dB 



= 0. 



x - d 2 n dXj 

^ dXidXj ~dB 



d 2 n 

dXidB 



which is an equation for the derivatives We assume 
this to be solved and thus the to be known. Inserting 
now A,- — A, 



B ■ into equation ^ we obtain 



n = n 



d 2 n ^ a 2 n dXj 

'dB 2 + ^ dXidB dB 

i 

from which the susceptibility is found as 
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(7) 



d 2 n 
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x - a 2 n dXj 

^ dX. t dB ~d~B 



In the presence of a magnetic field all single electron pa- 
rameters of the reference system have to be taken spin 
dependent, i.e. a hopping integral t — > (tf, tjj. It is then 

easy to see that the mixed second derivatives 



d 2 n 

dXii 



different from zero only for those A^ which are odd un- 
der sign change of the spin, that means quantities like 
(i| —£[.)■ At a nonmagnetic SP these are zero so the 
derivatives can be evaluated right at the SP. 
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